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The gravitational lens system MG J0414+0534 is formed by an elliptical galaxy at redshift z ~ 
0.96 and a quasar at z ~ 2.64. The system geometry is typical of lensing by an elliptical galaxy 
with the QSO close to and inside a fold caustic. It shows 4 images of the background source, 
and a partial Einstein ring is visible at optical wavelengths. It was observed with a global-VLBI 
array at A 18 cm in June 2008. We present here the imaging results and a preliminary lens model 
constrained by these observations. 
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1. The gravitational lensing system MG J0414+0534 

Predicted by general relativity, gravitational lensing (GL) occurs when light rays from a back- 
ground source travel through the gravitational potential of a foreground mass distribution. Multiple 
images of the same source are produced, with the observed image configuration (relative positions 
of the lensed images, fluxes ratios and time delays) carrying information on the mass distribution 
that is causing the effect. This phenomenon allows us to study mass distributions on different scales 
(corresponding to different regimes of lensing) and, since it is a purely gravitational effect, it has 
also been used as a cosmological tool to address questions related to dark matter halo properties, 
non-luminous satellite galaxies and the large-scale structure distribution in the Universe. 

The lens system presented here provides a beautiful example of an observed image config- 
uration due to GL by a foreground elliptical galaxy. It was discovered by the VLA-MIT survey 
(Bennett et al., 1986), and also later by other lens surveys (JVAS/CLASS) (Myers et al., 2003). It 
consists of a main elliptical galaxy lens (Schechter & Moore, 1993) at redshift z= 0.9584 ±0.0002 
(Tonry & Kochanek, 1999) and four images of a QSO at z= 2.639 ± 0.002 (Lawrence et al, 1995) 
separated by up to 2". Hubble Space Telescope (HST) observations (Falco et al., 1997) show a 
partial Einstein ring connecting the two brightest, merging images and a third image. The optical 
data also show the presence of a luminous satellite, 'X', more than 1" away from the main lens 
whose light distribution is not consistent with that of an elliptical galaxy (Schechter & Moore, 
1993). VLBI observations (Patnaik & Porcas, 1996; Porcas, 1998; Ros et al, 2000; Trotter et al., 
2000) have resolved the radio structures of the QSO images. All show a resolved core region and 
two jet-like structures (Ros et al., 2000; Trotter et al, 2000). The complex geometry seen with high 
angular resolution radio observations and the differences between the lensed images may confuse 
the lensing interpretation, but actually can be explained by magnification gradients and the source 
structure. On the other hand, given that the observed geometry depends on the properties of the lens 
mass, its complexity may be seen as providing detailed information on the mass distribution. This 
system also offers an interesting opportunity for CDM-substructure studies. None of the models 
cited above is able to reproduce the observed flux ratios between the lensed images, and Minezaki 
et al. (2009) infer a lower limit on the mass of a CDM subhalo which would explain the mis-match 
between the observed and modeled flux ratios, in particular the anomalous A1/A2 flux. 
The lensed quasar has also been found to host a powerful water maser (Impellizzeri et al., 2008). 
An accurate lensing model, constructed from high-sensitivity VLBI observations will be needed to 
interpret the unlensed properties of the water maser system. This was one of the motivations for 
the study presented here. 

2. New global VLBI observations at X 18 cm 

Utilizing all the information that the observed image configuration provides is not easy, how- 
ever, as the intrinsic (unlensed) source structure and the lens mass distribution need to be modeled 
at the same time. This can be done using the LensClean method (Wucknitz, 2004) which fits for 
the source structure in a non-parametric way and determines the best lens model based on it. The 
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results presented here provide additional detailed structural information on the lensed components, 
which is crucial in order to utilize the Lens Modeling method mentioned above. 

2.1 Observations and data reduction 

We observed the system at a wavelength of 18 cm on 7 June 2008 for 19 hours with a global 
VLBI array comprising the following 21 telescopes: [EVN] Effelsberg (100 m), Westerbork tied- 
array (equivalent to 93 m), Jodrell Bank (76 m), Onsala (25 m), Medicina (32 m), Noto (32 m), 
Torun (32 m), Urumqi (25 m), Shanghai (25 m), Hartebeesthoek (26 m), Arecibo 1 (305 m). 
[US NRAO] 9 of the 25 m dishes of the VLBA (Saint Croix, Hancock, North Liberty, Fort Davis, 
Los Alamos, Kitt Peak, Owens Valley, Brewster, Mauna Kea) and the VLA tied-array (equivalent 
to 130 m). Data were recorded in dual-polarization mode at 512 Mbps, distributed over 8 subbands 
per polarization, each of 8 MHz bandwidth. There were 16 frequency points (channels) per sub- 
band. The integration time was 1 second. The flux density of the target source is strong enough 
(Katz et al., 1997) so that phase-referencing was only needed to create an initial image prior to 
self-calibration. 

The data reduction and analysis were performed using the NRAO Astronomical Image Pro- 
cessing System (AIPS, version 31Decl0). The visibility phases were corrected for changes in 
parallactic angle, and a priori amplitude calibration was determined from measured system tem- 
peratures and antenna gains. Standard fringe-fitting was performed on the two phase calibrators 2 . 
Changes of gain as a function of frequency (across all the channels) were corrected via bandpass 
calibration. The spectral response of each telescope was analyzed from the data of the fringe- 
finders, and the final bandpass correction for each telescope was constructed using only the scan 
with the flattest amplitude and phase response. Subsequently, the two phase-reference sources were 
mapped by performing self-calibration (initially phase-only but finally phase and amplitude). The 
phase and amplitude corrections from the best models for both sources were applied to the target, 
which was then imaged. The data on the lens system were then phase self-calibrated with solution 
intervals of 4, 2 and 1 minutes. Throughout the whole process, the data were never averaged in 
frequency or time, in order to preserve a large field of view. This resulted in a rather large number 
of visibilities. 

2.2 Imaging results 

The J04 14+0534 lens system is 2" in size, and therefore wide-field techniques are required 
for imaging (Garrett et al., 1999). Instead of mapping this as a single field with milliarcsecond 
(mas) resolution (which would result in an image of quite unmanageable size) we used the AIPS 
task IMAGR, which provides a wide-field mapping procedure based on a multi-field approach. We 
mapped 3 sub-fields, centered on images B, C and the mid-point of images Al and A2. 

Our maps are shown in Fig. |l[ The central panel is a combination of optical and infrared HST 
images from the CfA-Arizona Space Telescope LEns Survey (CASTLES) (Falco et al., 1999) and 

'our target, at declination 5°, can be tracked for ~2 hours at Arecibo, close to the maximum of 2h:46m 
2 J0422 + 0219 and J0412 + 0438, 3.85 ° and 1.07 ° from the target; the second was used when Arecibo was in. 
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the outer panels are the new VLBI maps from the observations described here. A label in each 
panel identifies the lensed image. The restoring beam has a size of ~ 8 x 3 mas 2 and position angle 
(P. A.) ~ —9 deg. Previous high resolution VLBI observations of this system (Ros et al., 2000; 
Trotter et al., 2000) have shown the presence of a resolved core region and 2 jet-like structures in 
each of the 4 lensed images. Our new global VLBI data has higher sensitivity; a comparison of our 
maps with e.g. Fig. 2 of Trotter et al. (2000) shows that the structures are consistent with those 
shown in previous images and, moreover, in images Al, A2 and B new components in the outer 
part of the more elongated jet are clearly seen. We adopt the same notation as these authors to label 
the different components in the images: p and q for the two core components, r for the single jet 
component on one side of the core, s for the more extended jet regions on the other side and t the 
outer part of this. 

Image B shows a resolved core region extending over ~ 20 mas and two jet like features (r 
and s/t) not exactly aligned with it. This image is highly stretched in the direction tangential to that 
of the lensing galaxy. Images Al and A2 show the most complex structures, highly distorted in the 
tangential directions. Image Al has the two jet-like features parallel to each other and not aligned 
with the central core. The brighter jet (s/t) is highly elongated and extends for more than 80 mas in 
the NE direction. In image A2 the two extended jet-like features are also not aligned with the core. 
The brighter (s/t) is elongated in the NS direction over ~ 60 mas, while r is ~ 80 mas away in the 
SW direction. The core region extends over ~ 20 mas, showing a two-component structure. Image 
C shows a NS structure, with the core in between r and s/t; the northern jet component is brighter 
than the southern, and is more elongated. This image is the least bright (and hence least resolved) 
amongst the four. The total flux densities measured for each subfield agree with the flux ratios 
between the lensed images given by Katz et al. (1997). The lens models cited above predict that 
the two image pairs A1&A2 and B&C should each have opposite parities. The extended structures 
seen at radio wavelengths clearly confirm this although the lensing transformations are highly non- 
linear. 

3. Lens modeling 

The work done by Ros et al. (2000) to model this system, using X-band VLBI data as con- 
straints, has shown that the observed image configuration is reproduced quite well by modeling 
the main lens galaxy as a Singular Isothermal Ellipse plus a contribution from external shear, and 
assuming a Singular Isothermal Sphere for the luminous satellite X. 

Trotter et al. (2000) used C-band VLBA data to represent the gravitational potential as a Taylor 
expansion. They find that the mass distribution of the main lens galaxy is highly asymmetric, or 
a perturber is contributing to the lensing effect (see Section 5 in Trotter et al. 2000). Both studies 
find that the luminous satellite X needs to be included in the lens modeling, but they use a rather 
simple model for the lensed source. 

In Table [j] we show preliminary results from lens modeling which we have initiated, which 
will be used as a starting point for the more complex method mentioned above (Sec. [j]). As con- 
straints we have used the image positions of the centroid p and of feature r. 
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Figure 1: The gravitational lens system MG J0414+0534. Centre: CASTLES image. The 4 lensed images 
(Al, A2, B and C), the main lens galaxy and the luminous satellite X are visible. Outer panels: The radio 
VLBI maps from the new observations. In each panel the labels identify different regions of the source. 
The size of the restoring beam is 7.9 x 2.7 mas 2 in position angle -9.2 deg. The maps display contours 
of 1.4 mJy/beamx (-1,1,2,4,7, 16,32). The rms noise levels in the 4 images are 0.18, 0.18, 0.17, 0.09 
mJy/beam for Al, A2, B and C, respectively. 



Table 1: Best model parameters. The image positions of the core centroid, component r, and the optical 
positions of the main lens and object X were used as constraints. An accuracy of 0.1 mas for the relative po- 
sitions was assumed. The main lens galaxy was parametrized as a Singular Isothermal Ellipse with external 
shear and the satellite X as a Singular Isothermal Sphere. The lens strengths of the 2 lenses (Einstein radii 
b m /> bt), and the ellipticity and shear of the main lens (e, 7) were varied while searching for the best values 
for the ellipse and shear position angles (6 e , 9y). Afterwards all the parameters, together with the positions 
of the lenses, x m i, y,,,/ and x x , y x with respect to the core centroid in Al, were varied. 
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4. Summary and Outlook 

Our preliminary results from lens modeling confirm that the main lens galaxy and the luminous 
satellite X can reproduce the lensed image configuration for the system presented here. The values 
fitted for the ellipticity and sky position angle agree with those constrained by optical data (Falco 
et al., 1997). Analysis of these results, and error estimates on the model parameters are currently 
being done. 
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We plan to refine our results using a non-parametric lens-modeling approach following Wucknitz 
(2004); we will address the question concerning the radial mass distribution of the main lens galaxy, 
the position of the satellite X and any other dark substructure which is causing the flux mis-matches. 
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